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Abstract 

In this paper, we propose a model to investigate the electron thermionic emission from a single¬ 
layer graphene (ignoring the effects of substrate) and to explore its application as the emitter of 
thermionic energy convertor (TIC). An analytical formula has been derived, which is a function of 
temperature, work function and Fermi energy level. The formula is significantly different from the 
traditional Richardson-Dushman (RD) law for which it is independent of mass to account for the 
supply function of the electrons in the graphene behaving like massless Fermion quasiparticles. By 
comparing with a recent experiment [Kaili Jiang et al., Nano Research 7, 553 (2014)] measuring 
electron thermionic emission from a suspended single layer graphene, our model predicts that the 
intrinsic work function of a single-layer graphene is about 4.514 eV with a Fermi energy level of 
0.083 eV. For a given work function, a new scaling of T 3 is predicted, which is different from the 
traditional RD scaling of T 2 . If the work function of the graphene is lowered to 2.5 to 3 eV, and 
the Fermi energy level is increased to 0.8 to 0.9 eV, it is possible to design a graphene cathode 
based TIC operating at around 900 K or lower, as compared with the metal-based cathode TIC 
(operating at about 1500 K). With a graphene based cathode (work function = 4.514 eV) at 900 
K, and a metallic based anode (work function = 2.5 eV) like LaB6 at 425 K, the efficiency of our 
proposed-TIC is about 45%. 

PACS numbers: 79.70.+q, 79.20.Ws, 79.60.-i, 41.75.-i 


1 


INTRODUCTION 


Thermionic emission describes electrons evaporation from a heated cathode when the 
electrons gain sufficient energy from the thermal energy to overcome the potential barrier 
(or work function) near the cathode surface (see Fig. 1). The amount of current density 
J from the thermionic emission is determined by the Richardson-Dunshman (RD) law [Tj, 
given by 

J = AT 2 exp[-^\. (1) 

Here, T is the cathode temperature, $ is the work function of metal (independent of T), 
ks is the Boltzmann constant, the prefactor A = inernkg/h 3 = 1.2 x 10 6 A m -2 K -2 is the 
Richardson constant, e is the electron charge, and m is the electron mass. In general, RD 
law only works well for metallic-like materials, and a corresponding modification by using 
quantum models is required for wide band-gap materials and low electron affinity materials 

i- 

Since the mono-layer graphene [3] was exfoliated experimentally in 2004, many unique 
properties have been reported, such as linear band structure [4], ultra-high mobility (up to 
40000 cm-V _1 -s _1 ), and excellent conductivity [5]. Fundamentally, the linear band structure 
of graphene, the most intriguing property, makes it different from other 3-dimensional or 
bulk materials. For electron emission, it has been recently shown that the traditional emis¬ 
sion processes, such as field emission and photo-assisted over-barrier electron emission, may 
require further revisions (SHE] to account for the unique properties of graphene. As the crys¬ 
talline allotrope of carbon, the thermionic and field emission from carbon nanotube (CNT) 
had also been studied both experimentally [SHU! and theoretically [IMS!, which indicated 
that traditional emission models may not be valid for CNT. Recent experiments [T6l IT7 J have 
confirmed that RD law is not valid for thermionic emission from CNT. 

In this paper, we are interested to know if the RD law (Eq. [Tj) is valid for thermionic 
emission from a single-layer suspended graphene by assuming that the effect of the substrate 
is not important. Electrons in graphene behave as massless quasiparticles (m = 0), so the 
mass-dependent expression of the prefactor A in the RD law [Eq. (1)] is questionable as the 
supply function of the electron behaving like massless particles in graphene is not included. 
On the other hand, the electrons in graphene must exhibit a non-zero mass when they are 
collectively excited, which is only a few percents (0.01 to 0.03) of the intrinsic electron mass 
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Harvesting thermal energy is important to maintain a sustainable energy needs since 
nearly 60% of the energy input to our society is wasted as heat. One of the most common 
methods in converting heat into electricity is based on thermoelectric (TE) materials, which 
has progressed significantly since 1990s by using low-dimensional materials [19J . One of the 
limitation of TE power generators is low efficiency (< 40%) in the temperature ranges of 600 
to 1000 K (see Fig. 6 in Ref. pT9]) - Another method for the conversion of heat to electricity 
is known as the thermionic energy convertor (TIC), for which electrons are evaporated from 
a heated cathode into a vacuum and then condensed at a cooler anode [20]. Compared with 
TE based convertor (TEC), TIC normally operates at high temperature (above 1500 K) 
with higher efficiency (> 50%). However, it is difficult to operate TIC at low temperature, 
as common metallic cathodes having high work function cannot produce sufficient electron 
emission at low temperature. 

The other limiting factor of TIC is the space charge effects of the electrons traveling 
across the gap, which have been addressed either by reducing the gap spacing or by using 
Cs positive ions. Recently, some satisfactory progresses have been achieved to avoid this 
space charge problem by modifying the electric potential inside the gap m or by introducing 
heterostructures [22] . 

It is desirable if TIC can operate with relatively high efficiency (compared to TEC) at 
immediate temperature ranging from 700 K to 1000 K, which may significantly complement 
the ongoing research of TEC and also being alternative energy harvesting devices in this 
intermediate temperature range. Some recent approaches to improving the efficiency of TIC 
are achieved by using photo-enhanced thermionic emission from a wide-band gap emitter 
[23], engineering surface effects [24] and seeking alternative emitters [25] . 

Other than studying the validity of using RD law to describe the electron thermionic 
emission from graphene, we are also interested in exploring if a graphene-based cathode TIC 
can have better efficiency than a metal-based cathode TIC at similar operating condition 
like same work function. From our model, we first derive an analytic solution [Eq. (|8])] for the 
thermionic emission from graphene, which indicates that the traditional RD law has failed to 
describe thermionic emission from a single suspended layer of graphene [see Fig. 2(a)]. The 
calculated results are in a very good agreement with a very recent experimental data [26] 
[see Fig. 4(a)]. If the property of graphene can be tuned to some favorable parameters, we 
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also predict that the proposed graphene cathode based TIC can have higher current density 
(> 10 A/m 2 ) at 700 to 1000 K [see Fig. 3]. Finally, we show that a graphene based cathode 
(work function = 4.514 eV) at 900 K, and a metallic based anode (work function = 2.5 eV) 
like LaB 6 at 425 K, the efficieny of our proposed TIC is about 45 % [see Fig. 5]. 

THEORETICAL MODEL 

In Fig. la, we show the schematic diagram of our model. By assuming that the graphene 
based cathode has been heated up to a uniform temperature T c (without considering the 
heating process and uniformity issue for simplicity), the electrons will be emitted from the 
cathode to reach the anode biased at temperature T a , by overcoming the potential barrier 
as shown in Fig. lb. The work function of the cathode and of anode is, respectively, defined 
as 4> c and <f> a . The amount of the emitted current density from the cathode and anode is, 
respectively J c and J a , which will be calculated in this paper. 

The electron theory of graphene proposed by Wallace ra can provide the basis for 
practically all properties of graphene. This theory implies that electron in the graphene 
mimics the Dirac fermion, and that the equation of electron motion obeys the 2D Dirac-like 
equation. By using this model, the electron state is described by a two-component wave 
function for which the parallel electron energy is E p = hvf | k | in the low energy regime. 
The number of electron states per unit cell with energy between E p and E p + dE p is given 
by m 

p(E p )dE p = — J j d 2 k = -^-J^dE p , (2) 

where h is the reduced Planck constant and Vf (10 6 m/s) is the velocity of massless Dirac 
fermions in the graphene. The probability of an electron state with total energy E being 
occupied is given by the Fermi-Dirac (FD) distribution function, 

fF- D (E) = —- jj— - — w , ,, (3) 

1 + exp[(E - Ef)/k B T\ 

where Ef is the Fermi energy level (= 0.083 eV for intrinsic graphene). The number of 
electrons (per area and per time) perpendicularly crossing the graphene plane with total 
energy between E and E + dE and normal energy between E x and E x + dE x is given by 
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v x d 2 kdk x 


(4) 


j j 

= n fo 3 v 2 EpfFD(E)dE p dE x . 

In the derivation of Eq. (4), we have assumed that the normal energy component of 
electrons is E x = h 2 k x /2m. This assumption may be justified by the reasoning that graphene 
has atomic thickness for which the electrons are confined in a finite quantum well in the 
normal direction. The barrier height of the quantum well is assumed to be the intrinsic 
work function of a single layer graphene, which may be different from bulk graphite. Due to 
the quantum confinement, we solve the time-independent Schrodinger equation of electrons 
traveling in the normal direction to obtain the energy of the ground state, which is E x as 
given in the equations above. Using Eq. (4), we calculate the total number of electrons with 
normal energy between E x and E x + dE x at equilibrium condition, which gives 

j-rj roo 

N(E x )dE x = J e ( E ~ E x )f F . D (E)dE. (5) 

To solve Eq. ([5]) analytically, we make the follow assumptions, which will be verified by 
comparing the analytical result with numerical calculation as shown in Fig. 2(b). Firstly, 
we assume that only the electrons with energy greater than or equal to the work function of 
cathode (graphene) can be emitted, which means electron tunneling is completely omitted. 
This is reasonable as it is the definition of thermionic electron emission. It also implies that 
only high-energy component of the Fermi-Dirac distribution is important, for which we may 
replace the Fermi-Dirac distribution with the Maxwell-Boltzmann distribution /m- b(E) = 
exp[— (E — Ep)/kBT]. In doing so, Eq.(|5| is simplified to become 

N(E x )dE x = -JE-^exp[ ^-^-}dE x . (6) 

For thermionic emission, the current density of the emitted electrons along the direction 
perpendicular to graphene plane is calculated by 


pOO 

J(E f ,T)= / eN(E x )dE x . 


(7) 


In solving Eqs. (6) and (7), we obtain an analytical formula for thermionic electron emission 
from a single-layer graphene excluding the effect of substrate, which is 

ek%T 3 r $ — E f 


J(e f ,t) = 


V B J 

n h 3 v 2 


exp[- 


knT 


( 8 ) 
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It is clear that Eq. (8) is independent of mass m, solving the inconsistency of RD law 
having a finite mass term in the prefactor A, which is troublesome in applying RD law for 
electron emission from a single layer graphene. Note the work function is considered to be 
temperature-independent in our paper. Based on this interesting hireling and the importance 
of electron band structure determining the amount of electron thermionic emission, we 
speculate that all the materials with linear band structure may have the same scaling law 
reported here. 

It is common that an external and small voltage is used to collect thermionically emitted 
electrons at the anode. With the applied voltage, the potential barrier will be reduced to 
& e ff = — A<h (due to Schottky effect) as shown in the Fig. lb (dark blue line). Under 

this situation, the equation is modified as 

J(Ef,T) = p T 3 exp[- - ~ E ^~ A$ l, (9) 

where /3 = ek^/irffiv'j = 115.8 Am 2 K -3 , Ad* = (e 3 F/47re 0 ) 1 A represents the reduction in 
the potential barrier, and F (< 0.1 V/nm) is the electric held [29]. At higher held F > 1 
V/nm, quantum tunneling begins to dominate over the Schottky emission process, which 
requires new model [50], and it will not be addressed in this paper. Throughout this paper, 
we have assumed F — 0 unless it is specihed. 

To verify Eq.Q, it is important to have a more realistic surface potential prohle V(x) 
in order to estimate the contribution of electron tunneling through the potential barrier. In 
doing so, the classical image charge potential V irn is arbitrarily included in V{x) (near to 
the cathode), and V irn is [HI] 


lim ("^) 


e 2 . 1 nD 1 

47re D 2x ' (nD 2 ) — x 2 D 

n =1 v ' 


( 10 ) 


which can be simplihed to be [32] 


Um (•£) 


1.15a D 2 
x(x — D ) 


( 11 ) 


Here, a = e 2 ln2/87reo-D, D is the gap spacing, and e Q is the vacuum permittivity. The 
potential barrier V (x) in between the gap is 


V(x) = $ - V im (x) 


1.15 aD 2 
x(x — D ) 


( 12 ) 
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Note that the singularities at x = 0 and x — D can be avoided by restricting the positive 
values of x (within 0 < x < D) to have a positive V(x). This simplified image charge model 
had been verified by using a quantum image charge model [33], and it is valid for D > 10 
nm. 

To determine the effect of electron tunneling through the barrier V(x), we calculate the 
tunneling probability D(E X ) by solving 

D(E X ) = exp[-^-j^ J [(V(x) - E x )} 1/2 dx\ (13) 

« exp[-B{H - E x ) l '\ 

where H — $ — 1.15a/bln[ |^j 2 ] is the modified barrier height, b = a/ 1 — 4.6a/$ and 
B = 2(x 2 — xi)(2m) 1 / 2 /h. The values of x\ and x 2 are calculated by solving V(x) = 0. The 
amount of the tunneling current density is calculated by 

roo 

J T = e N(E x )D(E x )dE x (14) 

Jo 

_ 1 E j! 

~ Iwvj/e ^ ~ 2 ^ bT)3P ° lylog[3 ’ 

where polylog[3,x] is the polylogarithm function |34j . 

In the derivation of Eq. (14), we have assumed that only electrons around the Fermi 
level participate in the tunneling process. Based on the default parameters used in our 
calculation: Ef = 0.083 eV, = 4.514 eV and vj = 10 6 m/s, the term is nearly 

zero, which implies that the amount of the tunneling current density Jt can be ignored as 
compared with the thermionic current density given by Eq. ([8]). 

THERMIONIC EMISSION FROM GRAPHENE 

In Fig. 2(a), the thermionic emitted current density J based on Eq. ([8]) is plotted as a 
function of temperature T (log-log scale) for various $ = 3, 4, and 5 eV [solid lines] at fixed 
Ef = 0.083 eV, and Vf = 10 6 m/s. The results based on the RD law are also presented 
[dashed lines] for comparison. Here, we have assumed that the work function does not 
change with the temperature. For a given work function, a new scaling of T 3 is predicted 
by our model, which is different from the traditional RD scaling of T 2 . In the inset, we plot 
ln(J/T 3 ) (left y axis for our model) and ln(J/T 2 ) (right y axis for RD law) versus 1/T for 
work functions = 3 to 5 eV. 
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In Fig. 2(b), we compare the analytical formula Eq. (J8]) [red line] with the numerical 
solution of Eq. (Jt]) [circles] up to T — 2500 K (at Ef = 0.083 eV and $ = 4.514 eV). The 
comparison shows a good agreement up to T = 2000 K. The deviation at higher temperature 
is due to the approximation of Maxwell-Boltzmann distribution to Fermi-Dirac distribution 
in solving the integral. This comparison implies that Eq. ([8]) is accurate up to 2000 K (a 
practical temperature). 

In Fig. 3(a), the thermionic current density J [from Eq. (|8|] is calculated as a function 
of T (up to 2000 K) for various Fermi energy level Ef = 0.083 to 0.4 eV. From the figure, we 
see that increasing Ef can drastically enhance J. In order to decrease the temperature to T 
= 900 to 1000 K with J > 10 A/m 2 , it is required to have a higher Ef = 0.8 to 0.9 eV and 
also a lower work function $ = 2.5 or 3 eV as shown in Fig. 3b. Note that the calculated 
J is sensitive at around Ef = 0.8 and 0.9 eV and 4> = 2.5 to 3 eV. For example, if we use 
lower Ef = 0.8 eV and higher <3> = 3 eV, J is significantly to J < 1 A/m 2 (see the insert 
in Fig. 3b). It is worth to note that Fermi level can be tuned experimentally over a wide 
range (0.5 to 0.85 eV) via different methods such as chemical doping [35], and electrostatic 
gate voltage [36]. Some recent studies has also suggested that graphene can be engineered 
to have a lower work function (2.5 to 3 eV) [37]. 

There is a recent experiment measuring the thermionic electron emission from a single¬ 
layer graphene (suspended on a CNT him) in the range of T = 1600 to 1750 K [see Fig. 4 
in Ref. |26j]. In Fig. 4(a), we fit the experimental data [black square symbol] to our model 
with a work function of about $ = 4.514 eV and Ep = 0.083 eV, which shows a better 
agreement (black solid line) as compared with using traditional RD law with 4.66 eV for the 
bulk graphite (dashed blue line). In comparison, the RD law at 4.514 eV (dashed red line) 
shows much larger values and our model at 4.66 eV predicts much lower current values than 
the experimental data. This extracted work function (4.514 eV) is in excellent agreement 
with the experimental measured value around 4.5 eV [33], and also the theoretical value 
by DFT [39]. Theoretically, the reduction of work function is expected because when the 
thickness of bulk graphite is reduced down to single-layer graphene, the quantum size effect 
in the normal direction becomes very important, then the size-reduce-induced confinement 
causes the upward shift in electron total energy, and thus the work function of graphene 
become smaller [see Section 4.1 in the Ref. [30]]. 

On the other hand, if the graphene is placed on different metallic substrates, another 


recent experiment [25] has shown that the data is well described by the traditional RD law, 
as shown in the Fig.4(b). Here, the type of substrate will affect the amount of current 
significantly with a fitted and effective work function (based on RD law) ranging from 3.3 
to 4.6 eV including effects of the substrate [25]. The findings [reported in Figs. 4(a) and 
4(b)] indicate that the metallic based substrate will play an important role for thermionic 
emission from graphene. Note that the effects of the substrate have been completely ignored 
in our model, which is beyond the scope of this paper. 

THERMIONIC ENERGY CONVERTOR (TIC) 

By ignoring other energy losses, we calculate the ideal efficiency rj of a thermionic energy 
convertor (TIC) in using graphene as cathode (4.514 eV) at cathode temperature T c , and 
metallic material as anode of different work function $ a and temperature T a . The value of 
rj is calculated by 

T] 

Here, we have used Eq. ([8]) to calculate J c , and the RD law for J a . 

Figure 5 shows the calculated rj (solid lines) as a function <F a = 0.5 to 4.5 eV for different 
T a [K] = 300, 500, 800 and 1000 for a fixed graphene cathode property: T c = 1800 K, <3> c 
= 4.514 eV and Ef = 0.083 eV. From the figure, we see that, for a given T a , there is an 
optimal value of anode work function <3> a for which the TIC will have a maximal value of 
Vmax- As an example, we have rj max ~ 0.45 at T a = 800 K, T c = 1800 K, and <F a ~ 2.25 eV 
(blue line). This behavior remains valid for T c down to 900 K [see the inset]. The maximum 
efficiency r] max (left y axis) and its corresponding optimal values of anode work function $ a 
(right y axis) are plotted in the inset as a function of anode temperature T a = 300 to 800 K 
for two cathode temperature T c = 900 K (black dashed lines) and 1800 K (red solid lines). 
Using this insert, we can immediately know the maximum efficiency rjmax of TIC and obtain 
the required values of T a and 4> a . At low T c = 900 K (see black dashed lines in the inset), 
we will need T a ps 425 K and <l> a ~ 2.5 eV to have rjmax ~ 0.45. Note that the most recent 
work shows that LaB 6 heterostructure may have work function around to 2 to 3 eV [4Tj, 
and it may be a suitable candidate for the anode material for the proposed high efficiency 
TIC shown in Fig. 5. 


(Jc- Ja)($c-$a) 


Jc(&c + 2 ksTc) — J a (Q a + 2 ksT a ) 


(15) 
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CONCLUSIONS 


In summary, we have proposed to use single-layer graphene as the cathode material in 
the design of a high efficiency thermionic energy convertor (TIC). In order to provide some 
quantitative design parameters, we first develop a model to calculate the electron thermionic 
emission from a single-layer suspended graphene by ignoring the effect of the substrate. 
From our model, an analytical formula [Eq. ([8j)] is formulated, which has been verified with 
a numerical calculation and compared with experimental results. Our findings suggest that 
the traditional thermionic emission law governed by the well-known Richardson-Dunshman 
(RD) equation is no longer valid if the effect of the substrate is not important. Our model 
predicts a new scaling of T 3 , which is different from the classical RD scaling of T 2 . 

By considering that the graphene can be engineered to have a higher Fermi energy (around 
0.8 to 0.9 eV) [35j [36], and to have a lower work function (to 2.5 to 3 eV) [37], it is possible 
to have a graphene-based cathode to emit a sufficiently high current density (> 10 A/m 2 ) 
at T = 700 to 1000 K. 

A proposed TIC with a single layer graphene as cathode (work function = 4.514 eV) 
and metallic electrode as anode of different work function is designed for which the optimal 
anode work function and anode temperature are calculated to obtain a maximal efficiency. 
It is possible to have a TIC at an efficiency of about 45 percent with a cathode temperature 
of 900 K and with a metallic anode of about 2.5 eV work function (like LaB 6 anode) at 
around 425 K. These findings may shed new light on developing thermionic cathode using 
graphene and also thermal energy harvesting device like TIC. The work is also useful in the 
development of novel materials with low work function, like the new LaBe heterostructure 
with a lowering work function of 0.46 eV HD- 
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CAPTIONS 


FIG. 1: (Color online) (a) Schematic diagram of the electron thermionic emission process of a 
graphene-based thermionic energy converter (TIC), where T c and T a are the cathode temperature 
and anode temperature, respectively; J c and J a are current flow from cathode to anode and from 
anode to cathode, respectively, (b) Energy level for the TIC shown in (a). 4> c and 4> a are the 
work function of cathode and of anode, respectively. & e ff is the effective barrier height of the 
cathode due to Schottky lowering effect. The orange (and blue) line represents the potential profile 
without (with) Schottky lowering effect. The red line traces the thermionic emission process of an 
hot electron into vacuum. 

FIG. 2: (Color online) (a) The emitted current density J [A/m 2 ] as a function of temperature 
T [K] for our model (solid lines) and RD law (dashed lines) at three different work function 
[eV] =3, 4 and 5 (top to bottom). The inset shows ln(J/T s ) (left y axis) and ln(J/T 2 ) (right y 
axis) versus 1/T. (b) A comparison of the analytical formula (Eq. 8) [red line] with the numerical 
solution of Eq. (7) [circle] as a function of temperature (for Ef = 0.083 eV and 4> = 4.514 eV). 

FIG. 3: (Color online) (a) The current density J as a function of temperature T for various Ep 
[eV]= 0.083 (intrinsic), 0.1, 0.25 and 0.4 (bottom to top), (b )J as a function of T at higher Ep 
[eV] = 0.8 and 0.9, and lower 4> [eV] = 2.5 and 3. The inset presents the result at Ep = 0.8 eV 
and = 3 eV with a much lower J. 

FIG. 4: (Color online) (a) A comparison of our calculated results (black solid line: 4.514 eV, 
magenta solid line: 4.66 eV) and of RD law (red) [blue dashed line: 4.66 eV, red dashed line: 
4.514 eV] with experimental results [26] (black square) for electron emission from a single layer 
suspended graphene (no substrate effect). From the experiment, the emission area of graphene is 
5 x 2.5 nun 2 and the anode voltage (500 V) is separated from the graphene surface by 3.7 mm. 
The applied F is less 1 V / yrri so Schottky lowering effect can be ignored, (b) A comparison with 
experimental result [25] (symbols) and RD law (line), where the effect of substrate is important as 
the graphene is on top of a metallic substrate. 

FIG. 5: (Color online) Efficiency of a graphene-based cathode thermionic energy converter (TIC) 
as a function of the work function of metallic anode <h a at four different anode temperature T a 
= 300, 500, 800 and 1000 K and fixed cathotjlg temperature T c = 1800 K. The inset gives the 
maximum efficiency r] ma x (left) and the required 4> a (right) as a function of T a at T c = 900 K 



Model for thermionic emission 


(b) 



FIG. 6: (a) Schematic diagram of the electron thermionic emission process of a graphene-based 
thermionic energy converter (TIC), where T c and T a are the cathode temperature and anode 
temperature, respectively; J c and J a are current flow from cathode to anode and from anode to 
cathode, respectively, (b) Energy level for the TIC shown in (a). <3? c and <l> a are the work function 
of cathode and of anode, respectively. & e ff is the effective barrier height of the cathode due to 
Schottky lowering effect. The orange (and blue) line represents the potential profile without (with) 
Schottky lowering effect. The red line traces the thermionic emission process of an hot electron 
into vacuum. 
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FIG. 7: (a) The emitted current density J [A/m 2 ] as a function of temperature T [K] for our model 
(solid lines) and RD law (dashed lines) at three different work function <1? [eV] =3, 4 and 5 (top 
to bottom). The inset shows ln(J/T 3 ) (left y axis) and ln(J/T 2 ) (right y axis) versus 1/T. (b) 
A comparison of the analytical formula (Eq. 8) [red line] with the numerical solution of Eq. (7) 
[circle] as a function of temperature (for Ef = 0.083 eV and = 4.514 eV). 
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FIG. 8: The current density J as a function of temperature T for various Ep [eV]= 0.083 (intrinsic), 
0.1, 0.25 and 0.4 (bottom to top), (b ) J as a function of T at higher Ep [eV] = 0.8 and 0.9, and 
lower [eV] = 2.5 and 3. The inset presents the result at Ep = 0.8 eV and = 3 eV with a much 
lower J. 
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FIG. 9: (a) A comparison of our calculated results (black solid line: 4.514 eV, magenta solid 
line: 4.66 eV) and of RD law (red) [blue dashed line: 4.66 eV, red dashed line: 4.514 eV] with 
experimental results [26j (black square) for electron emission from a single layer suspended graphene 
(no substrate effect). From the experiment, the emission area of graphene is 5 x 2.5 mm 2 and the 
anode voltage (500 V) is separated from the graphene surface by 3.7 mm. The applied F is less 1 
V/gm so Schottky lowering effect can be ignored, (b) A comparison with experimental result [25] 
(symbols) and RD law (line), where the effect of substrate is important as the graphene is on top 
of a metallic substrate. 
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FIG. 10: Efficiency of a graphene-based cathode thermionic energy converter (TIC) as a function 
of the work function of metallic anode <h a at four different anode temperature T a = 300, 500, 800 
and 1000 K and fixed cathode temperature T c = 1800 K. The inset gives the maximum efficiency 
Vmax (left) and the required <3? a (right) as a function of T a at T c = 900 K (black) and 1800 K (red). 


20 
















